####################################################################################
###
###   Immigration, Public Housing and Support for the French National Front
###   Gloria Gennaro
###
###   Paper Table A2
###
####################################################################################


library(dplyr)
library(AER)


################################################################################
# Set Up
################################################################################

# Load Working Directories
source(paste0(wd_main, '/2_code/00_working_directories.R'))

# Load data and clean
source(paste0(wd_code, '/01_data_load_clean.R'))


################################################################################
# Estimation (Table A2)
################################################################################

modelli = list()

for (n in c(1,2,3)){
  
  print(paste0('-------------- Immigation tercile ', n, ' -------------- '))
  
  df1 = df_did[which(df_did$imm_quant_99==n),]
  temp = lm(data=df1, fn ~  m4 + m3 + m2 + p1 + p2 + p3 + running + 
               factor(CODGEO) + factor(policy_period) )
  r2 = summary(temp)$r.squared
  
  temp = coeftest(temp, vcov = vcovCL, cluster = ~df1$CODGEO)

  # Print output
    print(round(temp[1:6, 1:2], 3))
    print(paste0("Observations: ", attr(temp, 'nobs')))
    print(paste0("R-squared: ", round(r2, 3)))
    
  modelli[[n]] = temp
}


